A trapped surface in the higher-dimensional self-similar Vaidya spacetime

We investigate a trapped surface and naked singularity in a $D$-dimensional Vaidya spacetime with a self-similar mass function. A trapped surface is defined as a closed spacelike $(D-2)$-surface which has negative both null expansions. There is no trapped surface in the Minkowski spacetime. However, in a four-dimensional self-similar Vaidya spacetime, Bengtsson and Senovilla considered non-spherical trapped surfaces and showed that a trapped surface can penetrate into a flat region, if and only if the mass function rises fast enough [I. Bengtsson and J. M. M. Senovilla, Phys. Rev. D \textbf{79}, 024027 (2009).]. We apply this result to a $D$-dimensional spacetime motivated by the context of large extra dimensions or TeV-scale gravity. In this paper, similarly to Bengtsson and Senovilla's study, we match four types of $(D-2)$-surfaces and show that a trapped surface extended into the flat region can be constructed in the $D$-dimensional Vaidya spacetime, if the increasing rate of the mass function is greater than 0.4628. Moreover, we show that the maximum radius of the trapped surface constructed here approaches the Schwarzschild-Tangherlini radius in the large $D$ limit. Also, we show that there is no naked singularity, if the spacetime has the trapped surface constructed here.Comment: 13 pages, 5 figure

Self-intersecting marginally outer trapped surfaces

We have shown previously that a merger of marginally outer trapped surfaces (MOTSs) occurs in a binary black hole merger and that there is a continuous sequence of MOTSs which connects the initial two black holes to the final one. In this paper, we confirm this scenario numerically and we detail further improvements in the numerical methods for locating MOTSs. With these improvements, we confirm the merger scenario and demonstrate the existence of self-intersecting MOTSs formed in the immediate aftermath of the merger. These results will allow us to track physical quantities across the non-linear merger process and to potentially infer properties of the merger from gravitational wave observations

The runaway instability in general relativistic accretion disks

When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.Comment: 5 pages, 3 figures, minor corrections to match published version in MNRAS, +link to animatio

Accurate Evolutions of Orbiting Binary Black Holes

We present a detailed analysis of binary black hole evolutions in the last orbit and demonstrate consistent and convergent results for the trajectories of the individual bodies. The gauge choice can significantly affect the overall accuracy of the evolution. It is possible to reconcile certain gauge-dependent discrepancies by examining the convergence limit. We illustrate these results using an initial data set recently evolved by Brügmann et al. [Phys. Rev. Lett. 92, 211101 (2004)]. For our highest resolution and most accurate gauge, we estimate the duration of this data set's last orbit to be approximately 59MADM

New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions

We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators' spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions.Comment: 16 pages, 9 figures. The files with the coefficients for the derivative and dissipation operators can be accessed by downloading the source code for the document. The files are located in the "coeffs" subdirector

Numerical relativity with characteristic evolution, using six angular patches

The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular coordinate system has proved effective. Here we investigate the use of a six-patch system in characteristic numerical relativity, by comparing an existing two-patch implementation (using second-order finite differencing throughout) with a new six-patch implementation (using either second- or fourth-order finite differencing for the angular derivatives). We compare these different codes by monitoring the Einstein constraint equations, numerically evaluated independently from the evolution. We find that, compared to the (second-order) two-patch code at equivalent resolutions, the errors of the second-order six-patch code are smaller by a factor of about 2, and the errors of the fourth-order six-patch code are smaller by a factor of nearly 50.Comment: 12 pages, 5 figures, submitted to CQG (special NFNR issue

Dynamics and gravitational wave signature of collapsar formation

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formatio

Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.Comment: 7 pages, 16 figure

Dynamics and Gravitational Wave Signature of Collapsar Formation

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formation